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Abstract: In light of the recent solar and atmospheric neutrino data, we investi- 
gate the possibility of having an exactly degenerate spectrum for heavy right-handed 
Majorana neutrinos at the grand unification scale. The analysis is performed in the 
context of the minimal supersymmetric standard model with unbroken R parity and 
extended with three heavy Majorana neutrino fields in order to implement the seesaw 
mechanism. In the absence of a Dirac-type leptonic mixing, the only source of lep- 
j> ! ton flavour violation is the right-handed neutrino sector. Inspired by GUT-motivated 

relations among the quark, charged-lepton and Dirac neutrino Yukawa coupling ma- 
trices, and after the inclusion of the radiative effects, we determine the effective 
neutrino mass matrix at the electroweak scale. Using then the latest global analyses 
of the solar and atmospheric data at 99% C.L., we conclude that, within this frame- 
work, the only solar solutions compatible with the experimental data are the LOW 
and LMA solutions, being the latter the most favoured one. At 90% C.L., only the 
LMA solution is allowed. 
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1. Introduction 

The recent observation of solar and atmospheric neutrino deficits constitutes a solid 
evidence for physics beyond the Standard Model. The existence of such anomalies 
gives us a strong indication of neutrino oscillations among the three known families, 
which in turn imply that neutrinos have nonvanishing masses and nontrivial mixings. 
The latest atmospheric neutrino data from the Super-Kamiokande Collaboration [1] 
and CHOOZ experiments [2] tends to favour the — > v T oscillations with a large 
mixing angle # atm — 7r/4. On the other hand, the solar neutrino experiments [3]- [6] 
point towards the conversions v e — * z/ M and v e — > z/ r , and strongly favour the large 
mixing angle (LMA) MSW [7] solution over the small mixing angle (SMA), the LOW 
and the just-so vacuum oscillation (VO) solutions, as well as over the sterile neutrino 
hypothesis. 

From a theoretical viewpoint, modelling neutrino masses and mixings represents 
a challenge. First, the presently most favoured solution, namely the LMA MSW 
solution, is usually the most difficult one to achieve in realistic models. Second, 
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large mixing angles do not seem to fit naturally in our grand unification view. This 
is because in a very large class of string and grand unified theories (GUT's), the 
quark, charged lepton and Dirac neutrino mass matrices are typically related to each 
other and, therefore, one would expect small neutrino mixings in analogy to what is 
observed in the quark sector. The above reasoning assumes of course that the large 
neutrino mixing arises only from the Dirac neutrino mass matrix structure. Third, 
even if large neutrino mixings are consistently incorporated in the general picture of 
fermions, a natural question then arises, namely, why neutrino masses are so tiny 
compared to the rest of the fermions of the theory. 

A simple and economical way to overcome some of the above difficulties is to in- 
voke the so-called seesaw mechanism [8] which, in the presence of heavy right-handed 
Majorana neutrinos, u R , leads to an effective mass matrix for the light neutrinos. 
Rather than an assumption, the existence of the vrS is a natural requirement in 
the context of all GUT's with a group symmetry larger than SU(5). The seesaw 
mechanism by itself has, unfortunately, a very limited predictive power: the overall 
scale of the light neutrino masses is not uniquely determined, the mixing angles and 
the ratios of the effective neutrino masses are not fixed. 

Grand unified theories such as 5*0(10) [9] are a suitable framework not only to 
analyze fermion masses but also to implement the seesaw mechanism. One of the 
attractive features of the SO(10) model is that its gauge group is left-right symmetric 
and, consequently, there exists a complete quark-lepton symmetry in the spectrum. 
In particular, the fact that all left-handed (right-handed) fermions of each family fit 
into the single irreducible spinor representation 16 (16) of SO(10) and that the right- 
handed neutrino vr is precisely contained in this representation is remarkable. As it 
happens in all gauge groups, after the spontaneous symmetry breaking, the fermions 
will acquire masses out of the Yukawa couplings and the vacuum expectation values 
(VEV) of the Higgs fields. Several constraints on fermion masses are usually implied 
in these models [9]. For instance, if there is only one 10 Higgs multiplet responsible 
for the masses, then we have the relation Yu = Y fl = Y e = Y u (the indices 
U, D, £ and v stand for the up quarks, down quarks, charged leptons and Dirac 
neutrinos, respectively). Similarly, the existence of two 10 Higgs multiplets implies 
Yfl = and Yu = Y v . On the other hand, if the fermion masses are generated 
by a VEV of the 126 of S'O(IO), then the SU(4) symmetry yields the relations 
3 Yfl = — Y^ and 3Yu = — Y u . Of course, equalities such as the ones arising 
purely from the 10-dimensional representation cannot be exact in realistic models, 
since they imply undesirable relations among the quark and charged lepton masses. 
Additional assumptions are therefore necessary in order to predict the correct fermion 
spectrum [10, 11]. 

At this point it is worth recalling that the Majorana mass matrix, M R , for 
the right-handed neutrinos is usually much less constrained in unified models than 
the Dirac mass matrix. This has led many authors to discuss possible textures for 
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the heavy Majorana neutrino mass matrix within the seesaw framework [12]- [17]. 
A common feature to most of these analyses is the fact that the heavy Majorana 
spectrum turns out to be either hierarchical or partially degenerate. We may there- 
fore ask ourselves whether a complete degenerate spectrum is allowed for M# in this 
context. Similar studies have been previously performed in the case of degenerate 
left-handed Majorana neutrinos [18]- [20]. In particular, it has recently been argued 
[19] that if the Majorana masses of the light neutrinos are exactly degenerate at some 
high scale, then the one-loop radiative corrections are not sufficient to produce the 
required neutrino mass splittings and mixings at low scale. 

Although a fully degenerate mass pattern for the light neutrinos could be natural, 
if the Majorana masses arise from nonrenormalizable operators, and is welcome by 
cosmology, e.g. if relic neutrinos constitute the hot dark matter of the Universe, at 
a fundamental level it seems more plausible to construct and motivate patterns for 
the Dirac and right-handed Majorana matrices rather than for the effective neutrino 
mass matrix. In fact, the effective neutrino mass matrix is only a secondary output 
of the seesaw mechanism and, consequently, it may not be directly related to any 
symmetry of the theory. 

Of course, the observed hierarchy Am^ tm = |m| — m|| ^> Am^ = \m\ — m\\ 
among the light neutrino masses (i = 1,2,3) implies that one of the following 
mass patterns should be realized: |m lj2 | <C |m 3 | (hierarchical), |mi| ~ |m 2 | |m 3 | 
(inversely hierarchical) or \m\\ ~ |m 2 | ~ |m 3 | (almost degenerate). Thus, if one 
believes in the seesaw as the mechanism responsible for the smallness of neutrino 
masses, any realistic texture for the Dirac and Majorana matrices should lead to one 
of these patterns. 

In this paper we address the question whether right-handed neutrino degeneracy 
is compatible or not with the present experimental solar and atmospheric neutrino 
data. We assume complete degeneracy of the heavy Majorana neutrinos ur at GUT 
scale and also that the right-handed neutrino sector provides all the source of lepton 
flavour violation in the theory. Neglecting the rotation between the charged lepton 
and the Dirac neutrino mass matrices is a reasonable approximation, if the Higgs 
field has its components along suitable representations of 5*0(10) such as the ones 
mentioned above. This scenario is quite appealing since one starts from a situation 
where there is no leptonic mixing (similar to what happens in the quark sector where 
the mixing is described by the CKM matrix and is therefore small) to generate large 
mixing through the seesaw mechanism. 

When discussing neutrino masses and mixing angles, a crucial aspect is the run- 
ning of the various couplings from the unification scale down to the electroweak scale. 
In the presence of nonvanishing neutrino couplings, the renormalization group equa- 
tions (RGE) are modified. From the unification scale A x to the v R decoupling scale 
A R , the radiative corrections from the Dirac and right-handed Majorana neutrinos 
must be included. Below A R , when the seesaw mechanism becomes operative and 
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the vrs decouple from the theory, the running of the effective neutrino Yukawa cou- 
plings then becomes relevant. Radiative corrections could indeed provide potential 
contributions to the effective neutrino mass matrix or be even disastrous. A typical 
example is given by the RGE effects due to the charged-lepton Yukawa couplings 
[21]. 

The outline of the paper is as follows. In Section 2 we briefly present the general 
framework used in our analysis, with some attention given to the parametrization of 
the degenerate right-handed Majorana neutrino mass matrix. Section 3 is devoted to 
the study of the RGE effects from the GUT scale down to the electroweak scale. In 
order to account for the right-handed decoupling threshold at A#, the RGE analysis 
consists of two steps: first, the running of all the quark, charged-lepton as well as the 
Dirac and right-handed neutrino Yukawa couplings from Ax to A# and, second, the 
evolution of all the quark and charged-lepton couplings together with the effective 
neutrino mass matrix from Ar to the electroweak scale Mz- Such an analysis is 
required in order to verify whether the imposed patterns for the Yukawa coupling 
matrices at high scales are compatible or not with the low-energy data. In Section 4, 
the effective neutrino mass matrix, obtained at the electroweak scale, is confronted 
with the experimental data to determine which of the currently allowed solar neutrino 
solutions can indeed be realized in the present framework. Some numerical examples 
for the most favoured solution, namely the LMA MSW solution, and for different 
hierarchical patterns of the Dirac neutrino Yukawa couplings, are given in this section. 
Few comments on neutrinoless double (3 decay are also presented. These results allow 
us to motivate a simple analytical structure for the effective neutrino mass matrix in 
Section 5. Finally, we present our concluding remarks in Section 6. 



2. Theoretical framework 

We shall work in the context of the minimal supersymmetric standard model (MSSM) 
with unbroken R parity and extended with 3 heavy right-handed neutrinos. Super- 
symmetry (SUSY) is a natural choice to implement the seesaw mechanism. For 
instance, the hierarchy problem, usually present in a nonsupersymmetric seesaw, is 
automatically cured in the supersymmetric case. Indeed, the existence of very heavy 
right-handed neutrinos coupled to the Higgs field implies logarithmically-divergent 
radiative corrections to the Higgs mass ~ A 2 R YlH 2 ln(A^//z 2 ). While such correc- 
tions persist in nonsupersymmetric theories, they are naturally cancelled in SUSY 
by similar contributions coming from the right-handed sneutrinos. 

The relevant Yukawa mass terms are described by the Lagrangian (generation 
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indices are suppressed) 



C Y = QYuU c H 2 + QY D D c H ± + LY e E c ^ 

+ LY u v R H 2 + ^v R T M R v R + h.c. , (2.1) 

where Yjj , Yd , Y<? and Y u are the up-quark, down-quark, charged-lepton and Dirac 
neutrino Yukawa coupling matrices 1 , respectively; Mr is a 3x3 symmetric Majo- 
rana mass matrix which preserves the SM gauge symmetry; Q and L are the quark 
and charged-lepton left-handed doublets; U c , D c and E c are the up-quark, down- 
quark and charged-lepton right-handed singlets. The additional neutrino chiral fields 
v Ri (i — r) are right-handed Majorana fields responsible for the seesaw mecha- 
nism. Finally, Hi and H 2 are the hypercharge -1/2 and +1/2 MSSM Higgs doublets, 
whose neutral components acquire VEVs after the spontaneous symmetry breaking: 
(if?) = v cos , (H°) = v sin (3, with v ~ 174 GeV. 

In what follows we shall assume that Y„ oc and also that there is a left- 
handed alignment between Y^ and Y u at the unification scale Ax, which can be equal 
or higher than the decoupling scale A R , i.e. the scale at which the seesaw mechanism 
becomes operative and the v R s decouple from the light fields. As mentioned in 
the Introduction, these are reasonable assumptions in the context of grand unified 
theories. This means, in particular, that the only source of non-trivial leptonic mixing 
is the right-handed neutrino sector. In this scenario large neutrino mixings will be 
generated through the seesaw mechanism. 

Next we assume that the right-handed Majorana neutrinos are degenerate at 
Ax- In this case a simple and general parametrization [22] can be given for the 
right-handed Majorana neutrino mass matrix M R . Indeed, since M K is symmetric 
we can always write 

W T M R W = A R t , Y R = A^Mfl = W*W^ , (2.2) 

where W is a unitary matrix and Y R is a unitary symmetric matrix. It is then 
straightforward to show that the matrix Y R can be parametrized as 

Y R = He R<?> R-6» T , (2-3) 
where the rotation matrices He and are given by 

Re = I -s e c e , R^ = c 2(f> -s 24> . (2.4) 
\ 1/ V -s 24> -c 2 J 

Here we have introduced, for simplicity, the notations c v = cos ip , = sin <p . 
1 Wc assume all the Yukawa coupling matrices real. 
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Let us note that the parametrization (2.3) does not include the trivial case when 
CP is conserved and all the right-handed neutrinos have the same CP parity [22]. 
The CP- violating phase a turns out to be irrelevant in our further analysis and thus 
we shall set a = 0. In the latter case the complete form of Y R is given by 

/c 2 e + c 24> sj -suflsj -s e s 2(j> \ 
Y R = -s 2 esl s 2 e + c 2(j> c 2 e -c e s 2(j> . (2.5) 

\ -SeS 2(j) -CgS 2(f> -C24 ) 

Moreover, 

Mr- 1 = A^MRe^R/r 1 = A^R,R R/ = A^Y fl . (2.6) 

Notice that the parametrization (2.5) highly constrains the form of M R , which only 
depends on two angles 9, and the decoupling scale Ar . 



3. The effects of radiative corrections 

In this section we shall consider the evolution of the couplings from the unification 
scale Ax to the electroweak scale Mz- The relevant one-loop RGE, which include 
neutrino threshold effects, are summarized in the Appendix. Let us note that it is 
customary to take the SUSY breaking scale as A S usy — 1 TeV or A S usy — m t . Since 
in general the renormalization group effects that could possibly arise between Agusy 
and M z are negligible, we will not account for such effects in our analysis. Thus we 
set Asusy = M z . 

At the unification scale Ax the Yukawa matrices Y^ and Y u can be diagonalized 
by real and orthogonal matrices Ul,r and Vl,r such that 

U L ] Y e U R = diag ( y e (A x ) , y M (A x ) , vA^x) ) , 

Vi t Y,VR = diag(y 1 (Ax),y 2 (Ax),l/3(Ax)) , (3.1) 

where y^(Ax) denotes the respective Yukawa couplings at Ax- The hypothesis of 
left-handed alignment between Y u and Y^ implies Ul — Vl- In this case, Y v and Yi 
can be simultaneously diagonalized, leaving the charged-current Lagrangian invari- 
ant. This means that there is no "Dirac-type" leptonic mixing at Ax- Obviously, 
this diagonalization process does not affect the spectrum of the right-handed Ma- 
jorana matrix and, therefore, the resulting matrix M.' R = Vl M R V L can also be 
parametrized in the form of Eq. (2.5) with a different set of the parameters 9 and (f>. 

3.1 Running Y v from Ax to A R 

Our first step in discussing the effects of radiative corrections on neutrino masses and 
mixing angles will be the study of the evolution of the Dirac Yukawa coupling matrix 
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Y v from GUT scale to the right-handed neutrino decoupling scale. After performing 
the diagonalization given in Eq. (3.1) we have 



6! 0' 

Y„(A A ) =. ( y,(A. v ) \ e 2 

1 



(3.2) 



where the parameters 



2/i (Ax) _ y 2 (A x ) 

e l — 71 7 i e 2 



1/3 (Ax) ' 1/3 (Ax) 



(3.3) 



have been introduced. 

The running of Y u from Ax to A R is governed by Eq. (A. 4). From this equation 
and neglecting all but the third-generation Yukawa couplings 2 we get 



1 

Y„( V„) / | 1 | Y„(A.,v) . 
,00/,. 



(3.4) 



where 
Ji = exp 



16 7T 2 



T u -G u )dt 



tx 



exp 



16 7T 2 



3y t 2 + vl- - 3#2 ) * 



(3.5) 



/„ ~ exp 



16 7T 2 



tR 



3y 2 3 + y 2 T )dt 



tx 



(3.6) 



with t x = \n(A x /M z ) and t R = \n(A R /M z ) . 
From Eqs. (3.2) and (3.4) we find 

/cf 0\ 

Y„(A fl ) = K D , A D = y 3 (Ax) /i 4 , < 2 = J" 1 e 1>2 . (3.7) 
\0 1/ 

3.2 Running M# from Ax to Ar 

The running of the right-handed Majorana mass matrix from Ax to A R is de- 
scribed by Eq. (A. 10) given in the Appendix. This equation allows us to relate the 
right-handed neutrino mass matrices M,r(A,r) and M^(Ax): 



\M R (A R )]ij = [Mfl(Ax)]ij exp 



1 



8vr 2 



tx 



(y 2 + y 2 ) dt 



(3.8) 



2 This is a reasonable assumption considering that we will focus on cases where yi <C 1/2 <C 2/3- 
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In the limit where only the third Dirac neutrino Yukawa coupling y 3 is kept, we find 



1 \ (\ 

M R (A R ) =(010 \I,,.(.V V ) 1 | (3.9) 
OOIrJ \00 I r , 



where 



I R ~ exp 



_8tt 2 

The inverse of M R at A R is then given by 



y 3 dt 

tx 



(3.10) 



no o ^ 




fio \ 


10 


|Y.| 


10 


\0 0IrJ 




Voo/^V 



(3.11) 



where Y R has the form given in Eq. (2.5). 
3.3 M^, ff from high to low scales 

In order to study the low-energy consequences of our assumptions, imposed at high 
scales, we must find the effective neutrino mass matrix M£ ff at the typical neutrino 
experiment scale, i.e. the electroweak scale (~ Mz). If R parity is an exact symme- 
try in the MSSM, the lowest (5-dimensional) operator in the superpotential which 
violates lepton number and generates a Majorana mass for the left-handed neutrinos, 
v, is given by [23, 24] 

W vv = ^kv t v Hi Hi + h.c. , (3.12) 

where k is a 3x3 symmetric matrix. The effective neutrino mass matrix is obviously 
given by 

Mf = k(H° 2 ) 2 = Kv 2 sin 2 (3 . (3.13) 

The supersymmetric seesaw mechanism is realized by adding to the superpoten- 
tial the following terms: 

W = W^mssm - \v c r Mr v c r + u R T Y u L H 2 . (3.14) 

Below A R one can consider the effective superpotential given by 

W cS = W MS sM + l(Y l/ LH 2 ) T M R 1 (Y l/ LH 2 ) , (3.15) 
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which is obtained by integrating out the heavy right-handed neutrino fields vr in 
Eq. (3.14). The second term in W e g will originate in the Lagrangian a Majorana 
mass term of the type 

1 



C vv = — v 1 kvH 2 u H 2 " + h.c. , 
2 



where 



K 



(3.16) 



(3.17) 



Let us note that at this stage the matrices k and M£ ff are defined at the de- 
coupling scale Ar. To compare them with the experimental ones, obtained at the 
electroweak scale, we must run these matrices down to Mz through the corresponding 
RGE. Taking into account Eqs. (3.7) and (3.11) we have 



K(A i? ) = A^ 1 Y^(A i? )Y i? Y,(A i? ) 



K 



D 



A 



R 



e{ 
e' 2 
I R l 



R 



e[ 
e' 2 
I R l 



(3.18) 



The running of the matrix k from Ar to Mz is described by the RGE given 
in Eq. (A. 14). From the latter equation, and neglecting the charged-lepton Yukawa 
couplings and y e , we get the following relation between k(Mz) and k(Ar): 



k{M z ) 



1 
1 
I r 



«(Afl) 



1 
1 
h 



where 
h = 



exp 



16 7T 2 



tz 



T K — G K ) dt 



tR 



~ exp 



16 7T 2 



tz 



tR 



(3.19) 



6 

6y t 2 - -g\ - 6gi I dt 



(3.20) 



exp 



16 7T 2 



Using Eq. (3.18) we finally obtain 



with 



k{M z ) 



K e 



5 1 0^ 



A 



R 



1 



tz 



y 2 T dt 



tR 



'8i 0' 



5 2 Yr(A x ) 5 2 



1 



$1,2 — hff Ir* e l,2 , hff-dRd,, 1 , 



(3.21) 



(3.22) 



(3.23) 



K en = yl{A x )lll 2 I 2 T l;? n 



(3.24) 
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and l\ , I2 given in Eqs. (3.5) and (3.20), respectively. 

From Eqs. (2.5) and (3.22) we obtain our final expression for the neutrino effective 
operator at the electroweak scale: 



k(M z ) = -5! 5 2 s 2e sl 81 (sj + c 2<t> c 2 e ) -82 c e s 2(f) . (3.25) 

\ -81 s e S24> -8 2 C0S24, ~ c 24> I 

Let us remark that the mass spectrum and mixings arising from this matrix are 
invariant under the transformations (9 —> , — > it — (p) and (0 — > 7r — , — > 0), 
and, obviously, under any transformation resulting from their successive application. 
This simply corresponds to a phase redefinition of the left-handed neutrino fields. 
Thus, from now on we restrict ourselves to the intervals < 0, <fi < tt/2. 

In the limit where the evolution of the third Dirac neutrino Yukawa coupling y 3 
and the r Yukawa coupling y T can be neglected, we find the approximate expressions: 

• >Htt • <-> 

and, consequently, 

At this point few remarks are in order. First, as it has already been noticed 
by several authors [21], the parameter I T ~ 1. This can be easily checked by using 
either the exact expression (3.21) or the approximate expression given in Eq. (3.26). 
Moreover, the parameter I e ff, which controls the RGE effects of the Dirac neutrino 
Yukawa coupling ratios (cf. Eq. (3.23)), is also close to the unity. This can be seen 
from Fig. 1 where we present the dependence of I e ff on the value of the scale \x 
for the two initial conditions Y u = Y{/ and Y u = 3Yj/. We note that I e ff ~ 1 
for a wide range of tan/5. Notice also that for the case Yj, = 3Y[/, the use of the 
approximate expression for I e ff, obtained by neglecting the running of the Yukawa 
couplings, is not well justified. This is because in the latter case y 3 (A x ) = 3y t (A x ) 
and the evolution of the couplings becomes more relevant. 

4. Confronting the neutrino mass matrix with the experimen- 
tal data 

The analysis of the neutrino oscillation data indicates that, contrary to what happens 
in the quark sector, where the mixing is always small, the mixing angle is large 
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H (GeV) 



Figure 1: The evolution of the RGE coefficient I e ff = IrI^ 1 from Ax to A# , assuming 
two different Dirac neutrino textures, Y# = Yj/ and Yc = 3Yj/ , and for the values of 
tan/3 = 10 and tan/3 = 55. The solid curves correspond to the exact expressions as defined 
by Eqs. (3.6) and (3.10). The dashed lines correspond to the approximate curves as given 
in Eqs. (3.26). 

in the 2-3 neutrino sector. The latest atmospheric neutrino data from the Super- 
Kamiokande Collaboration [1] and CHOOZ experiments [2] implies a large mixing 
angle ^atm — 7r/4. In what concerns the 1-2 neutrino sector, there are four different 
scenarios depending on the solar neutrino problem solution. 

Since the present experimental data constrains the neutrino mass squared dif- 
ferences and the mixing angles, it is necessary to diagonalize the effective neutrino 
mass matrix Mf at the electroweak scale. The neutrino mixing matrix U, which 
relates the flavour and mass eigenstates, 



such that U T M^ ff (M^) U = diag(mi ,m 2 , m 3 ), where m ; are the light neutrino 
masses; c^- = cos 0y , = sin 9ij . 

We also remark that the lepton mixing matrix can have, besides a Dirac-type 
phase (analogous to that of the quark sector), two other physical phases associated 




(4.1) 



can be parametrized in the standard form [25] 




(4.2) 
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Atmospheric and reactor neutrinos 





Aml tm (eV 2 ) 


tan 2 6> atm 


\U e3 \ 






(1.1 - 7.3)xl(T 3 


0.3 - 3.8 


< 0.27 






[1.4-6.1]xlCT 3 


[0.4 - 3.1] 


[< 0.23] 






3.1 x 10 3 


1.4 


0.07 




Solar neutrinos 




Ami ( ey2 ) 


tan 2 6> 


r = Am|/Am 2 tm 


LMA 


(2.0 - 20)xl0- 5 


0.1 - 0.7 


(0.03- 1.8)xl0" 


-l 




[2.3 - 10]xl0~ 5 


[0.2 - 0.5] 


[0.04 -0.7] x 10" 


l 




4.2 x 10 5 


0.26 


0.014 




SMA 


(0.4- l)xl0" 5 


(0.2 - 2)xl0~ 3 


(0.55-9.1)xl0 _ 


-3 


LOW 


(0.5-2)xl0~ 7 


0.5-0.9 


(0.07- 1.8) xlO" 


-4 


Just-So 2 


(5-8)xl0~ 12 


0.5-2 


( 0.68 -7.3) x 10" 


-9 



Table 1: Constraints on neutrino masses and mixing angles at 99% C.L. [90% C.L.] 
coming from global analyses of the solar, atmospheric and reactor neutrino data [27, 28]. 
The numbers in bold correspond to the best-fit values. 



with the Majorana character of neutrinos. If CP is conserved these phases are 
either or it. In the analysis that follows we shall assume all these phases vanishing. 
However, as we shall see in Section 4.2, such phases play a significant role in processes 
such as the neutrinoless double (3 decay [26] . 

4.1 Neutrino oscillation data 

The relevant parameter set for the study of solar and atmospheric data is usually 
identified with the quantities 

Am 2 & = Am 2 2 = m\ - m\ , Am 2 tm = Am 2 3 = m 2 - m\ , 

tan 2 # = tan 2 9 12 , tan 2 atm = tan 2 6 23 , U e3 = sin 9i 3 , (4.3) 

where all mixing angles are assumed to lie in the range [0, 7r/2]. In Table 1 we summa- 
rize the present constraints on neutrino mass squared differences and on the mixing 
angles obtained from the recent global analyses [27, 28] of the solar, atmospheric and 
reactor neutrino data at 99% C.L. [90% C.L.] . 

Apart from the overall factor K e ff/A R , the effective neutrino operator k(M z ) 
defined in Eq.(3.25) depends on the angles 9 and 0, which parametrize the heavy 
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right-handed neutrinos, and the parameters Si : 2 . For the structure of k to be compat- 
ible with the low-energy data, we should take into account not only all the constraints 
on r = Am|/Am^ tm and the mixing angles, but also the experimentally allowed 
region for the up-quark mass ratios, since within our assumptions, these ratios fix 
the parameters 5\ and 5 2 at A^: 

The above ratios can be determined by running the RGE from the electroweak 
scale Mz to the unification scale Ax, which we assume to be Ax ~ 2xl0 16 . Since 
these ratios are quite insensitive to the value of tan f3 in a wide region, we shall take 
tan/3 = 10 as the input value in our analysis. Using the experimentally allowed 
values for the quark masses at the electroweak scale [29, 25], the following ranges are 
then obtained at A X - 

2.8xl(T 6 < 5 1 < 1.5xl0~ 5 , 
9xl(T 4 < 5 2 <4xl0" 3 . (4.5) 

In our numerical study we randomly varied all the parameters in the ranges: 
< 5i,2 < 1 , < 9 < vr/2 , < (f) < tt/2 . In Fig. 2 we plot the allowed region 
in the plane (S u 5 2 ) for the three different MSW solar solutions (LMA, SMA, LOW) 
at 99% C.L. . We assume that the right-handed Majorana neutrinos are exactly 
degenerate at the unification scale Ax- The dashed area delimited by the dotted 
lines corresponds to the allowed region for the Dirac neutrino coupling ratios, as 
defined by Eqs. (4.5) and assuming Y u oc Yjj at Ax- The intersection of the filled 
areas with the dashed area provides us with the points (^i,^) compatible with the 
low-energy data and the initial conditions imposed at Ax- From the plots it is clear 
that the SMA solution is excluded. The LOW solution is barely compatible with the 
data at 99% C.L., and certainly excluded at 90% C.L. . Finally, the LMA solution 
is the most favoured one, and at 90% C.L. it is the only one that survives. 

The results in the parameter space (0, 9) are shown in Fig. 2d at 99% C.L. . We 
notice that the allowed region for LMA is more extensive than the SMA and LOW 
regions, and the SMA region is in turn wider than the LOW one. This is a direct 
consequence of the hierarchy r Q LMA > r SMA > r LOW . For the VO (Just-So 2 ) solar 
solutions an extreme fine-tuning is required between 9 and in order to get a solution 
compatible with the solar and atmospheric neutrino data. The latter solutions are 
therefore not presented in our plots. 

Since the LMA solution is the most favoured in the present framework, the rest 
of our discussion will be devoted to it. In Figs. 3 we show the allowed region at 99% 
C.L. for the LMA MSW solution, in different planes of the input parameters S 1: 5 2 , 4> 
and 9. A remarkable feature is the fact that for small values of Si and 5 2 , which are 
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required in our framework, one has ~ 7r/4 and 9 ~ n/2. We remark that the same 
conclusion for the value of can be drawn by looking at the structure of the effective 
neutrino mass matrix (3.25), since for small values of 5i >2 , the relation cos 20 <C 1 
is required to avoid a quasi-diagonal pattern in the mass matrix. 

Let us conclude this section by giving few numerical examples (see Table 2). 
We consider two different GUT-motivated relations for the Dirac neutrino Yukawa 
coupling matrix: Y u = Yjj and Y v = SYjj. In both cases, consistency with the 
present LMA solar, atmospheric and reactor neutrino data was achieved. It was not 
possible however to obtain the best-fit values, simultaneously for all the parameters. 
By slightly modifying the initial assumption yi/2/3 = Uu/ut , we were then able to 
reproduce the best-fit values for all the solar and atmospheric neutrino data, as shown 
in the last column of the table. Notice also that the obtained value \U e3 \ ~ 0.05 is 
close to the best-fit value \U e3 \ = 0.07 given in Table I. 

As suggested by our full numerical analysis (see Figs. 3), the angles 9 and were 
chosen to be close to rc/2 and 7r/4, respectively. The degeneracy imposed at Ax for 
the right-handed Majorana neutrinos, \M\\ = \M 2 \ = \M 3 \ , is lifted in all cases by 
radiative corrections. At the decoupling scale A^ we have \Mi\ ~ |M 2 | ~ Ir^-r < 
|M 3 | ~ Ar. For the effective neutrinos with masses TOj , we obtain at the electroweak 
scale a hierarchical spectrum with 



|toi| < \m 2 \ c± \J Am\ 2 < \m 3 \ ~ y Am^ . (4.6) 

Finally, we also notice the stability of the atmospheric and solar neutrino parameters, 
r v = Aml 2 / Aml 3 , tan 2 9i 2 , tan 2 9 23 and U e3 with the energy scale. In fact, from 
A R to M z no significant running of these parameters is observed. This is related 
with the fact that I T ~ 1. 

4.2 Neutrinoless double (3 decay 

The Majorana nature of massive neutrinos can be traced by looking for processes 
where the total lepton charge L changes by two units, i.e. AL = 2. The information 
coming from these type of events can be crucial not only to establish the nature of 
neutrinos but also to reveal some aspects about CP violation in the leptonic sector, 
to which neutrino oscillations are insensitive. In particular, if neutrinos have a Dirac 
signature, then the two extra Majorana phases which appear in the neutrino mixing 
matrix will be absent. Moreover, even if CP invariance holds, some information 
about the relative CP parities and neutrino mass spectrum can be extracted. One 
example of such processes is the neutrinoless double (3 decay ((3(3 0u - decay) of the 
type 

(A,Z) -> (A,Z + 2) + 2e~ , (4.7) 
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Figure 2: Allowed region in the plane (81,82) for the three different MSW solar solutions 
(LMA, SMA, LOW) at 99% C.L., if we assume Y u oc Yjj and an exactly degenerate 
spectrum for the right-handed Majorana neutrinos at the unification scale Ax- The dashed 
area corresponds to the allowed region for the up quark mass ratios at Ax- At 90% C.L., 
the only solution that is compatible with the experimental data is the LMA one. The 
corresponding allowed region in the plane (<fi,6) is shown in Fig. 2d at 99% C.L. . 



which can occur through the exchange of massive Majorana neutrinos coupled to 
electrons. The decay rate is proportional to the so-called "effective Majorana mass 
parameter" 

|<m>| = \mi U^+m2 U^ 2 + m 3 U^ \ . (4.8) 
Although no experimental evidence has been found so far for this kind of processes, 
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Figure 3: Allowed region at 99% C.L. for the most favoured solar solution, namely the 
LMA MSW solution, in different planes of the input parameters 61,82, (ft and 6, within our 
framework. For small values of the parameters <5i and 82, the angles (ft and approach ir/4 
and 7r/2, respectively. 



an upper bound has been obtained by the 76 Ge Heidelberg- Moscow experiment [30]: 

I < m > I < 0.35 eV [ 90% C.L. ] . (4.9) 
If CP is conserved, then Eq. (4.8) can be written as 

I < m > I =| rfi P mi jj -2 + v cp m2 ^2 + ^cp mg ^2 | ^ m . > , (4.10) 

where r]f p = ±1 are the neutrino CP parities. By examining the invariants of the 
effective neutrino mass matrix given in Eq. (3.25), one can conclude that one of the 
neutrinos should have a negative CP parity and the other two should have positive 
CP parity. This is verified for the numerical examples considered in Table 2, where 
the neutrino with mass m 2 has ^ p = — 1, while the other two have rjf^ = +1. In 
this case, Eq. (4.10) reads 

|<m>| = \miU£ -m 2 U^ + m 3 U^\ . (4.11) 

A particular feature in the above expression is the possibility of a cancellation among 
the different terms [31, 26], such that |<m>| ~ 0. In fact, for all the numerical 
cases considered in Table 2, such a cancellation indeed occurs. This can be easily 
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0.34 


0.26 
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0.47 


0.57 


1.4 


Pe3\ 


0.18 


0.18 


0.05 


At M z = 91.2 GeV 


Keff 


0.31 


2.18 


0.30 


mi (eV) 
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1.9xl0~ 3 


1.6X10" 3 


m 2 (eV) 


(-) l.OxlO- 2 
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m 3 (eV) 
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3.9xl0~ 2 


5.6xl0~ 2 


Am 2 2 (eV 2 ) 
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1.2xl0~ 4 


4.3xl0~ 5 


H ( e V 2 ) 


1.5X10" 3 


1.4xl0~ 3 


3.1xl0~ 3 


r„ 


6.7xl0~ 2 


8.1xl0~ 2 


1.4xl0~ 2 


tan 2 0i2 


0.34 


0.34 


0.26 


tan 2 023 


0.47 


0.57 


1.4 




0.18 


0.18 


0.05 


<m>| (cV) 


2.0xl0- 12 


2.4xl0~ 12 


9.8xl0- 15 


m u (MeV) 


1.9 


2.0 


2.2 


m c (GeV) 


0.68 


0.69 


0.66 


m t (GeV) 


174 


174 


175 


Down-quark masses 


rn d = 4.61 MeV 


m s = 91.8 MeV 


m b = 2.99 GeV 


Chargcd-lepton masses 


m e = 0.49 MeV 


m M = 102.8 MeV 


m T = 1.75 GeV 



Table 2: Numerical results for tan/3 = 10. We assume right-handed neutrino degeneracy 
at Ax, which is then lifted by radiative corrections. Two different Ansatze were considered 
for the Dirac neutrino Yukawa coupling matrix, Y v = Yu and Y v = SYu- By slightly 
modifying the initial condition yi/ys = y u /vt it is possible to reproduce the best-fit values 
for all the solar and atmospheric neutrino data, as shown in the last column of the table. 
The negative CP parities are indicated by (— ). 

understood if we recall that the effective Majorana mass parameter is nothing but 
the absolute value of (Mf) n . Therefore, from Eq. (3.25) we have 

|<m>| oc | 8\ (c 2 9 + c 2(j> s 2 e ) | . (4.12) 

Since 6 ~ rr/2 and ~ rr/4, in order to account for the low-energy neutrino data, 
and 5i is small, it is now clear from Eq. (4.12) why |<m>| is highly suppressed in 
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our framework. 



5. A simple analytical form for Mf 

In this section we shall present a simple analytical example which can reproduce the 
results obtained in the last section. We have seen that, in the context of complete 
degeneracy of right-handed neutrinos and imposing proportionality between Y v and 
Y{/, the only possible solution at low energies is the LMA. Moreover, as it is clear 
from Figs. 3 this would require 9 ~ 7r/2 and ~ 7r/4. On the other hand, if one 
assumes in Eq. (3.25), = 7r/4 and 9 = n/2, then 



K (M Z ) = -pL 5\ . (5.1) 




In this case the neutrino masses are 



Keff 



{m ll m 2 ,m 3 } = f i ^ sin 2 (3 { 6 1 , 5, , S 2 2 } . (5.2) 

Taking into account Eqs. (5.1) and (5.2), it is straightforward to show that, inde- 
pendently of the value of the parameters 5± and 62, one always has sin 2 (2 6> 2 3) = 0, 
which is not phenomenologically viable. 

At first sight the above result may seem contradictory, since on one hand our 
numerical study clearly shows that values of 9 ~ 7r/2 and <fi ~ 7r/4 are indeed 
compatible with the low-energy neutrino data, but on the other hand, a simple 
analytical check gives a completely opposite answer. As we shall see shortly, the 
reason for this mismatching lies in the fact that small deviations, such as 59 = 
n/2 — 9 <^ 1 and 5(p = <p — 7r/4<^l, and/or small perturbations from our original 
assumptions, turn out to be crucial. To show that this is indeed the case, let us relax 
the condition of proportionality between Y u and for the smallest Dirac neutrino 
Yukawa coupling by putting 

e 1= ae , e 2 = e = — — — , (5.3) 

and write 9 and as 

<P = l + be\ 0=|-e, (5.4) 
where a and b are some positive real constants. From Eq. (3.23) we have then 
5i = aleffl' 1 e 2 ~ ae 2 , 5 2 = I e ff -C 1 e — e • 
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The effective neutrino mass operator (3.25) will be approximately given in this 

case by 



/a 2 e 4 (l-26) -ae 2 -a 
K (M z )^^lle 2 \ -ae 2 l-26e 4 -l). (5.6) 

Ar \ -a -1 26 

Due to the smallness of the parameter e we have, to a very good approximation, 

K(M z )~^e 2 ^ , s'=| 1 -1 | . (5.7) 

The invariants of the matrix k! are 

T = A! + A 2 + A 3 = 1 + 26, 

X = Ai A 2 + Ai A 3 + A 2 A 3 = -1 - a 2 - 2 6 , 

A = AiA 2 A 3 = -a 2 , (5.8) 




where Aj are the eigenvalues, 



A 1 = 2 v /pcos(!)+T, 



with 



A 2 ,3 = -2VPcos(-±-j+T, (5.9) 



2T 3 Tv T 2 v a 

q = H - A , p= ^, cos« = q — . (5.10) 

H 27 3 '^93' p 3 / 2 V ; 

The corresponding eigenvectors are given by 

Vj = — [ a ( 1 — Aj ) , — Aj , Aj ( Aj — 1 ) ] T , (5.11) 

Til 

where Hi are the normalization coefficients. 

In Fig. 4 we present the region in the (a, 6)-plane, which is compatible with the 
low-energy solar and atmospheric neutrino data at 90% and 99% C.L. for the LMA 
solution (cf. Table 1). We see that no fine-tuning is required in order to reproduce 
the data. To obtain the best-fit point (black dot in the figure), we assume at the 
unification scale Ax the same values for the charm and top Yukawa couplings as in 
the last column of Table 2, i.e. y c (Ax) = 1.7xl0~ 3 and y t (Ax) = 0.64 . Taking then 
as input parameters 

y c (A x ) _ CKwin - 3 



2.65xl0~ 3 , a = 0.15, 6 = 0.34, = 6.3x10 GeV , (5.12) 



19 



1.2 




0.4- 



0.2- 




90 % C.L. - 



99 % C.L. 



0.0 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



a 



Figure 4: Allowed region for the parameters a and b of the triangular texture (5.7) for 
k(Mz), in the case of the LMA MSW solution and considering the neutrino data at 90% 
and 99% C.L. summarized in Table 1. The black dot corresponds to the point which leads 
to the best-fit values of the solar and atmospheric data as shown in Eqs. (5.12)-(5.13). 

the following values are obtained at the electroweak scale for tan (3 — 10 : 



which correspond to the best-fit values in the neutrino mixing plane for the atmo- 
spheric neutrinos [27] and the LMA solar neutrino solution [28]. 

6. Conclusions 

In this paper we have considered the possibility of having an exactly degenerate 
spectrum for heavy right-handed Majorana neutrinos at the grand unification scale. 
Our theoretical framework was based on the minimal supersymmetric standard model 
with unbroken R parity and extended with three heavy Majorana neutrino fields. 



K eff = 0.3 , 
m 1 = 1.6xl0 -3 eV , 
Am 2 12 = 4.2xl0~ 5 eV 2 , 
tan 2 12 = 0.26 , 



<?i = l.OxlO" 6 , 5 2 = 2.7xl0- 3 , 

m 2 = 6.7x 10~ 3 eV , m 3 = 5.6x 10~ 2 eV , 

Am 2 3 = 3.1xl0- 3 eV 2 , r v = 1.4x 10~ 2 , 

tan 2 # 23 = l-4, |L7 e3 | = 0.05, 



<m>\ = 9.3xl0" 15 eV 



(5.13) 
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We assumed no Dirac-type leptonic mixing at GUT scale, and thus the right-handed 
neutrino sector constitutes the only source of lepton flavour violation in this context. 
We have then studied the renormalization group effects, which include the neutrino 
threshold effects, from high to low energies. 

Taking into account the allowed range for the quark mass ratios at GUT scale 
and the effective neutrino mass matrix obtained through the seesaw mechanism, we 
have investigated which of the possible solutions to the solar neutrino problem can 
be accommodated in our framework. Inspired by GUT-motivated relations among 
the quark, charged-lepton and Dirac neutrino mass matrices, and based on the latest 
global analyses of the solar and atmospheric data at 99% C.L., we have then con- 
cluded that the only solar solutions compatible with the experimental data are the 
LOW and LMA solutions. In fact, at 90% C.L. only the latter is allowed. Even in 
this case, the obtained values for the neutrino oscillation parameters are far from the 
best-fit points. In this sense, right-handed neutrino degeneracy embedded in GUT 
scenarios seems to be disfavoured by the low-energy neutrino data. We emphasize 
however that, by slightly relaxing the imposed relations among the Yukawa coupling 
matrices at high scale, one can easily reach the best-fit values. Our assumption of 
degeneracy is in fact a very restrictive one, since the right-handed Majorana neutrino 
matrix is parametrized in terms of only two angles and an overall scale A#. On 
the other hand, a fully or partially hierarchical mass spectrum for right-handed neu- 
trinos can be easily achieved with a much less constrained structure, and one could 
therefore expect that such patterns are easier to be reconciled with the low-energy 
neutrino data. 

Finally, we are also aware of possible contributions to lepton flavour violation 
processes like /i — > ey . As it has been recently pointed out by Casas and Ibarra 
[32], in supersymmetric theories with seesaw- like neutrino masses the prediction for 
the branching ratio of this process is, in general, larger than the experimental upper 
bound. In particular, scenarios such as the one considered here, where top-neutrino 
unification occurs, are severely constrained unless the leptonic Yukawa coupling ma- 
trices satisfy very specific requirements. 
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A. Renormalization group equations 



In this Appendix we present the relevant one-loop renormalization group equations 
for the MSSM extended with three right-handed Majorana neutrinos. 

Between Ax and A R the evolution of the up-quark, Yu , down-quark, Y D , 
charged-lepton, Yp , and Dirac neutrino, Y v , Yukawa coupling matrices is given 
by 



Ifi 2<*Yt/ 

16 7T = 

dt 




-G u ) + (3Y u Y\ 7 + Y d Y] 3 )' 


Yu, 


(A.l) 


o dY d 
16n 2 

dt 


~(T D 


- G D ) + (3Y D Y^ + YuY 




Y D , 


(A.2) 


dYp 
16tt 2 , = 
dt 


~(Te- 


-G e ) + (3Y e Yl + Y u Yi)~ 


Y e , 


(A.3) 


i6 y Yu = 

dt 


\t v 


-G v ) + (3Y v Yl + Y e Y\) 


Y 




(A.4) 



where t = ln(/i/Mz), \i is the renormalization scale. The quantities Tj and G{ are 
defined as 



Tu 


= Tr(3Y l/ Y t [/ + Y,Yt) 


, Gu 


= y£ 3 + 3 fe + 


13 2 
15 ^ ' 


(A.5) 


T D 


= Tr(3Y D Y^ + Y,Y{) 


, Gd 


16 2 . o 2 . 

= y 9 3 + 3 ^2 + 


7 2 
15*' 


(A.6) 


T e 


= Tr(3Y D Yl 3 + Y,Yj) , 


G e = 


92 + 15 9l ' 




(A.7) 


T v 


= Tr(3Y l/ Yj / + Y,Yt) . 


G v ~- 


= 3 g 2 + - g x , 




(A.8) 



where g$, g 2 , and g\ are the SU(3), SU{2) and C/(l) gauge couplings, which satisfy 
the RGE 



lQ7T 2 ^£ = b ig f , h = (33/5, 1, -3) . 
The evolution of the right-handed neutrino mass matrix, Mr , is given by 



(A.9) 



Sir 2 



,dM 



R 



dt 



M R (Y u Ylf + (Y u Yl) M R . 



(A.10) 



Below Ar , the renormalization group equations for the up-quark and charged- 
lepton Yukawa couplings get modified since after decoupling of the right-handed 
neutrino states, the contribution of Y„ should vanish. In this case the RGE for Y^ 
and Yi read 



16 7T 



,dY, 



v 



16 7T 2 



dt 
,dYe 

dt 



T U -G U ) + (3Y U Y\ J + Y D Y^ D ) 



T £ -G,) + 3Y,Y{ 



Ye, 



(A.H) 
(A.12) 
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with 



T f ; = 3Tr(Y [7 Y^). (A. 13) 

Finally, the evolution of the effective neutrino mass matrix n from to Mz is 
given by 

16 vr 2 ^ = ( T K - G K )k + ( Y,Yj ) k + k ( Y £ Yj ) t , (A.14) 

where 

T K = 6Tr(YuYl) , G K = Qg 2 2 + \gi. (A.15) 

o 
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